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We investigate color superconducting phase at high density in the extended Nambu-Jona-Lasinio 
model for the two flavor quarks. Because of the non-renormalizability of the model, physical observ- 
ables may depend on the regularization procedure, that is why we apply two types of regularization, 
the cut-off and the dimensional one to evaluate the phase structure, the equation of state and the 
relationship between the mass and the radius of a dense star. To obtain the phase structure we 
evaluate the minimum of the effective potential at finite temperature and chemical potential. The 
stress tensor is calculated to derive the equation of state. Solving the Tolman-Oppenheimer-Volkoff 
equation, we show the relationship between the mass and the radius of a dense star. The dependence 
on the regularization is found not to be small, interestingly, dimensional regularization predicts color 
superconductivity phase at rather large values of /i (in agreement with perturbative QCD in contrast 
to the cut-off regularization), in the larger temperature interval, the existence of heavier and larger 
quark stars. 

PACS numbers: ll.lO.Kk, ll.30.Qc, 12.39.-x 



I. INTRODUCTION 

Variety of new phases of the matter consisting of 
quarks and gluons is expected in the high density QCD. 
One of them is the broken color S'[/c(3) symmetry phase 
where the color superconductivity takes place, Refs. [1-4] 
Quark-quark interaction is attractive in the color anti- 
symmetric 3 channel. This force destabilizes the filed 
Fermi sea leading to the Cooper instability. Then a com- 
posite operator constructed from a color anti-triplet di- 
quark pair develops a non-vanishing expectation value. 
For the case of the two flavor quarks the diquark con- 
densation breaks the color SUc{^) symmetry down to 
the 5[/c(2). This broken symmetry phase is called two- 
flavor color superconducting (2SC) phase. Because of the 
high degeneracy of states near the Fermi surface, such 
a symmetry breakdown behavior is a subject of a non- 
perturbative description in QCD. 

To investigate non-perturbative QCD effects we are 
working in the low energy effective theory, namely in the 
Nambu-Jona-Lasinio (NJL) model in which the chiral 
symmetry is broken down dynamically The model 
is useful to evaluate low-energy phenomena in hadronic 
phase of QCD. The NJL model is extended in a simple 
way to include the attraction in the 3 qq channel. It is 
based on the point-like four-fcrmion interactions between 
quarks. Since the four-fcrmion interaction is the dimen- 
sion 6 operator, the model is non-renormalizablc in four 
space-time dimensions, therefore some regularization is 
needed to obtain finite results. 



Most of analysis have been using momentum cut-off 
regularization, where the cut-off scale is determined phe- 
nomenologically. Unfortunately this cut-off scale often 
breaks some symmetries of the model. Moreover the crit- 
ical chemical potential where the color superconductivity 
takes place is of the order of the cut-off scale (ultraviolet 
cut-off may even hit the Fermi sea cut-off). In such a sit- 
uation, it is expected that the regularization procedure 
has a non-negligible effect on the analysis of the color su- 
perconductivity. That is why in the present paper we an- 
alyze the extended NJL model by using the dimensional 
regularization as well. In the dimensional regularization 
the space-time dimension is analytically continued to less 
than four. 

Some high density states are observed in astrophysical 
objects. In the core of dense stars, like neutron stars in- 
terior, quark stars and so on, the color superconducting 
phase is expected to take place. Much attention has been 
paid to such dense stars to find an evidence of the color 
superconductivity [1, 01. Characteristics of the dense 
stars such as the mass and the radius are of great in- 
terest. A constraint on the mass and the radius of the 
star can be found by solving the Tolman-Oppenheimer- 
Volkoff (TOV) equation This solution depends on 
the equation of state (EoS), i.e. on the relationship be- 
tween the energy density and the pressure, of the quark 
matter inside the stars. It was pointed out that the ex- 
istence of the color superconducting phase may decrease 
the minimum of the radius of such dense stars d, 01 , [9- 
29]. The minimal size of the stars is analyzed precisely. 
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A finite contribution from the color neutrality and the (3 
equilibrium are analyzed in [11-15]. Nucleon and vector 
meson contributions are calculated in [16-19]. Rotating 
compact stars are considered in . 

In the paper the extended NJL model is regarded as 
a low energy effective theory of QCD. We consider two 
flavor quarks and investigate properties of compact stars 
in the 2SC phase. In Sec. II we introduce the model 
Lagrangian and the fermion propagator in color super- 
conducting phase. In Sec. Ill the phase structure of the 
model is evaluated at finite temperature and chemical po- 
tential. We numerically calculate expectation values of 
the composite operators (i/'^A) and (ip'^tp) and the quark 
number susceptibility. In Sec. IV we evaluate the energy- 
momentum tensor and determine the EoS . In Sec. V we 
solve the TOV equation with the account of the obtained 
EoS and show the relationship between the radius and the 
mass of a dense star in the color superconducting phase. 
Finally we give concluding remarks. 



II. EXTENDED NJL MODEL 

In a state with large chemical potential the 3 diquark 
channel interaction plays essential role in creation of a 
color Cooper pair of quarks. To study color supercon- 
ductivity one can extend the NJL model to include the 3 
diquark interactions explicitly. The temperature, T, and 
the chemical potential, /i, is introduced to the theory via 
the imaginary time formalism. The model is defined by 
the Lagrangian density 

+ Gs{{^^pf + {i^aoilhTjktpakf} (1) 
+ GD{i^''aj£jk£ldl5'llJdk){i'ipfieirne-'}gl5'llj%rn), 

where the indices a,b,d,e, f,g and j,k,l,m denote the 
colors(l,2,3) and flavors (m, d) of the quarks, m is the 
quark mass matrix, diag(mu,md), Tjk represents the 
isospin Pauli matrices, V'^ is the charge conjugate of the 
fleld V'- Gs and Gd are the effective coupling constants 
for the qq scalar and the diquark channel respectively. 
The third line in Eq. ([T]) exposes the attractive force in 
the 3 channel of the quark-quark interaction. 

For practical calculations it is convenient to introduce 
the auxiliary fields: scalar cr, pseudo-scalar tt and di- 
quark A*", to write down the Lagrangian density as 
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where the Nambu-Gor'kov representation [35| is used. 



and G represents the quark propagator, 

= 

/ — m — i^74 — (T — Z757r • t — zA°£e°75 
\ —iA^ ££^75 — m + ifij4 — cr — i757r • 



(3) 
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From the equation of motion for the auxiliary fields we 
get the following correspondence. 
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The generating functional for the Lagrangian density 
^ is given by 

Z = y"I?(*,§,A^A^^a,7^)e^'*'"(^— 
= M y"I?(A^A''^(T,7r)detl/2G-l 
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where A/" is a constant and the sources are 
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The quark mass affects the scalar a channel [361] , we set 
TV ~ Q. SUc{i) symmetry allows one to set A^'^ = and 
A = A'^. Here we study the phase structure for the case 
of the two flavor quarks. 



III. PHASE STRUCTURE 

The phase structure is determined by the effective po- 
tential. We set Jif = Jvj> = in the generating functional 
dHl) to obtain the effective potential. 
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where denotes the volume of the system, ^ = J d^~^Q 
Due to the space-time translational invariance the ex- 
pectation values of a and A are independent of the 
space-time coordinates. Assuming the isospin symmetry, 
niu — TTid = m and fiu = fJ-d = fJ-, the effective potential 
reads 
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where trl= 



2^/2, /3 = l/T and 



^ = v/p2 + (a 



^± = E±fi, E± = \/e± 
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(12) 
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Since the third term of the right hand side of Eq. (fTTj) 
involves a divergent integral, we use the dimensional reg- 
ularization (2 < D < 4) to regularize it [sl, [111. After 
integration over the angle variables, we have 



K=ff(a, A) 
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with A = 4-\/7r/(27r)^/^. For the cut-off regularization, 
the third term of the right hand side is obtained by the 
replacement 
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where A is a cut-off scale. The first three terms of the 
integral in Eq. (|13p are naively divergent. The divergent 
part of the Eq. p3p in the dimensional regularization 
reads 



dp p^~^ {E+ +E-+E) 



dp , 



r(-f)r(^) 
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The following renormalization conditions are used: 
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where Mq is a renormalization scale. 

The parameters of the model are fixed to reproduce 
some observables at T = /i = 0. Using the pion 



mass rriTT = 138 MeV, the pion decay constant /tt = 
92.4 MeV and the quark mass m = 4.5 MeV, we ob- 
tain A = 720 MeV and Gs = 3.80 x 10"^ MeV"^ for 
the cut-off regularization. Here we use the relationship, 
Gd — (3/4) Gs, corresponding to the one gluon exchange 
interaction. In the chiral limit, the critical temperature 
becomes 170 MeV at /i = with these parameters. 

In the dimensional regularization one can not avoid 
renormalization which introduces an additional parame- 
ter, renormalization scale, Mg. In Ref. [s^l the scale is 
simply fixed at the dynamically generated fermion mass. 
The similar procedure is also adopted in [3§|. Here we 
fix all parameters by reproducing physical observables. 
Reproduction of the critical temperature places an ad- 
ditional constraint on the parameters. We fix the pa- 
rameters for the dimensional regularization to reproduce 
the pion mass, the pion decay constant, the quark mass 
and the critical temperature 170 MeV at m and 
^ = 0. Then we have D = 2.28, Mq = 127 MeV, 
G's = 0.603 MeV^"^ and G^ = (3/4)G^ for the dimen- 
sional regularization. In Ref. [40| the renormalization 
scale is fixed by using the decay width of the 7r° 27. 
It gives the renormalization scale near the constituent 
quark mass. 

If the auxiliary field A = A'^ develops a non-vanishing 
value, the diquark condensation takes place. It implies 
that the QCD gauge symmetry is broken down from 
SUc{?>) to SUc{2) and the 2SC phase is realized. The ex- 
pectation values of a and A are found by observing the 
minimum of the effective potential. To find the minimum 
we numerically evaluate the effective potential . The 
expectation values are shown as functions of the chemi- 
cal potential in Figs. 1 and 2 near a typical temperature 
of dense star, T — 1 keV, in the cut-off and the dimen- 
sional regularizations respectively. We observe similar 
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FIG. 1: Behavior of the expectation values of auxiliary fields 
a and A for m = and m = 4.5 MeV at T = 1 keV in the 
cut-oiT regularization. 

behavior of (cr) and (A) in Figs. 1 and 2. As is known, 
the current quark mass enhances the chiral symmetry 
breaking and contributes to the crossover behavior, i.e. 
the expectation value, (cr), smoothly disappears as the 
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FIG. 2: Behavior of the expectation values of auxiliary fields 
a and A for m = and m = 4.5 MeV at T = 1 keV in the 
dimensional regularization. 



chemical potential fi increases. At the massless limit, 
the critical chemical potential for the chiral symmetry 
breaking is found to be /i ~ 300 and 220 MeV in the case 
of the cut-off and the dimensional regularizations respec- 
tively. A large mass gap is observed in Fig. 2. There is 
no coexistence phase with non-vanishing (a) and (A) in 
the dimensional regularization. The current quark mass 
hardly affects the behavior of the expectation value (A) . 
The behavior of (A) in Fig. 2. is quite different from that 
in Fig. 1. As is clearly seen in Fig. 2, the absolute value 
I (A) I monotonically increases as a function of fi. In the 
cut-off regularization, (A) is suppressed near the cut-off 
scale, ^ ~ A (= 720 MeV). 

We also calculate the quark number susceptibility Xq 
m mi m m, m which is defined by 



FIG. 3: Quark number susceptibility for /i — 11, 17, 23, 25.5, 
29 MeV in the model with cut-off regularization. 



X,^-|^n((a),(A)), 
where n((cr), (A)) is the quark number density, 
n((a),(A)) = (^V>=-|-Kff((a),(A)) 
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Figures 3 and 4 show the behavior of Xq as a function of 
the chemical potential in the cut-off and the dimensional 
regularizations, respectively. In Fig. 3 we show the be- 
havior of the susceptibility Xq near the phase boundary. 
A coexistence phase can be realized and different critical 
points are observed for (a) and (A) in the cut-off regular- 
ization. In Fig. 3 the maximum of at /i ~ 312 MeV for 
T = 29 MeV is induced by the chiral symmetry break- 
ing. Decreasing the temperature, moves this peak to the 
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FIG. 4: Quark number susceptibility for ^ - 
65.8, 67.8, 70.8 MeV in the model with dimensional reg 
ularization. 



right until T - 23 MeV. The gap ai ^jl ~ 327 MeV for 
T = 29 MeV is generated by the color symmetry break- 
ing. This gap moves to the left as the temperature de- 
creases. The susceptibility for the peak (gap) higher than 
that for the gap (peak), implies a crossover (the first or- 
der phase transition). The chemical potentials for the 
peak and the gap degenerate at T ~ 17 MeV, it corre- 
sponds to the critical end point (CEP). There is a simpler 
behavior of Xq in the dimensional regularization. There 
is no coexistence phase as is seen in Fig. 4. Such a picture 
is consistent with the Fig. 2. 

The phase structure of the chiral S'C/l(2) x SUr{2) 
symmetry and 2SC in the cut-off and the dimensional 
regularizations are shown in the Figs. 5 and 6, respec- 
tively. We plot the boundary of the phase where Xq has 
the peak or the gap for m = 4.5 MeV and also the bound- 
ary where (a) and/or (A) disappear at the massless limit, 
TO — 0. The coexistence phase induces a complex struc- 
ture around the CEP in the cut-off regularization. We 
find the tricritical point (TCP) at (T,^) = (46.7,281) 
and CEP at (T,^) = (17.0,314). As is clearly seen in 
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Fig. 5, the 2SC is suppressed as /i approaches to the cut- 
off scale. In Fig. 6 we observe TCP at (T, fi) = (87.2, 194) 
and CEP at (67.8, 226). In the dimensional regulariza- 
tion the 2SC gets enhanced as ii increases. 
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FIG. 5: Phase diagram for the model with the cut-off regu- 
larization. The solid, the dashed and the dashed-dotted lines 
denote the first, the second order phase transition and the 
crossover, respectively. We also plot the secondary maximum 
of Xq by the dotted fine. 
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FIG. 6: Phase diagram for the model with the dimensional 
regularization. The solid, the dashed and the dashed-dotted 
lines denote the first, the second order phase transition and 
the crossover, respectively. 



IV. ENERGY-MOMENTUM TENSOR 

Here we discuss the EoS in the extended NJL model. 
The energy-density and the pressure of the system are 
given by the energy-momentum tensor. Thus we calcu- 
late the energy-momentum tensor at finite T and fi in the 
dimensional regularization and compare the result with 
that in the cut-off regularization. The energy-momentum 



tensor is defined by 
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where the parenthesis in the subfix indicates symmetriza- 
tion. 

Following the imaginary time formalism, we introduce 
the temperature and the chemical potential. We perform 
Wick rotation and substitute the Lagrangian density ([2]) 
into (I20l). 
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To study phenomena in the static star, we take the 
thermal average of the energy-momentum tensor. It is 
given by the expectation value in the ground state where 
the equations of motion. 
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are satisfied. Therefore the expectation values of the 
stress tensor elements are simplified to 

(r44) = Yiu,-tv[-fi{idl-iti){^{x)^v)) 
1 ' -(A)^ 
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To regularize the composite operator with the quark 
fields at the same point we consider quark fields at the 

thtly different space-time points, x and y [H, S^, 
The expectation values of the composite operators 
{'tp{x)ip{y)) and {ij}'^ {x)ip'^ [y)) are the diagonal matrix el- 
ements of the quark propagator in the Nambu-Gor'kov 
representation. 

Substituting the fermion propagator obtained in Sec. II 
into Eqs. ([^5)1 and and integrating over the angle 
variables, we obtain 
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where Mq is the renormalization scale defined in the con- 
ditions and (fTT]) . It should be noted that the stress 
tensor is renormalized to have the correct mass dimen- 
sion in the four dimensional space-time. Taking the four 
dimensional limit and applying the modification (|14p , we 



numerically evaluate the integral in Eqs. (|27p and ([28)1 in 
the cut-off regularization. 

Since the divergent parts of the Eqs. (|27p and (I^S)) do 
not depend on T and fi explicitly, we can evaluate the 
divergent parts in the limit, fi and /3 — > oo. Using 
the dimensional regularization, we evaluate analytically 
the momentum integral. Because of the general covari- 
ance in this limit, one derives the same value for (T44) 
and (Tu), 
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where (cr) and (A) depend on T and 

Thus we obtain a finite expression for the energy- 
momentum tensor 
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We plot the energy density, — (T44), in Fig. [7] and the 
pressure, {Tu), in Fig. [51 The origin of the energy density 
and pressure can not be fixed in our model. Here we set 
(r44)|T=;.=o = and {Tu)\t=^=o = by hand. In 2SC 
phase, (A) > 0, — (T44) and {Tu) monotonically increase 
as functions of /i. We find the equation of state from 
these results. It plays an essential role in the study of 
the structure of dense stars. As is shown in Fig. 9, there 
is a difference between the dimensional and the cut-off 
regularizations. A harder state is observed in the cut-off 
regularization for smaller —{T44). On the other hand, a 
harder state is obtained in the dimensional regularization 
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FIG. 7: Behavior of the stress tensor —{Ta) as functions of 
the chemical potential at T = 1 keV. 



FIG. 9: Equation of state, i.e., relationship between —{T44) 
and {T,,). 
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FIG. 8: Behavior of the stress tensor (Tu) as functions of the 
chemical potential at T = 1 keV. 



for larger —{T44). It affects the core structure of dense 
stars with 2SC phase inside. 

Next we examine the thermodynamic relationships. In 
the thermodynamics the energy density p is derived from 
the generating functional, Z, through the following dif- 
ferentiation, 



(32) 



The pressure of the system, P, is obtained by differenti- 
ating Z with respect to $7, 
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In the present approximation for the generating func- 
tional the energy density and the pressure become 
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First we consider the expression for the energy density. 
Substituting the explicit expression for the effective po- 
tential (fT3|) into the first term of the right hand side of 
the Eq. dSll), we get 
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The second term of the right hand side of the Eq. (|34p is 
the number density (|19p . times the chemical potential /x. 
Therefore the energy density p is given by 
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This expression exactly reproduces —{T44) of the 
Eq. ((27l) . Thus the time component of the stress ten- 
sor, — satisfies the thermodynamic relationship for 
the energy density ([M]) . 

For the space component of the stress tensor, {Tu), the 
situation is not so simple. After partial integration overp 
in the effective potential Eq. (HH]) or (Tu) of the Eq. ([28l). 
we have 
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At the limit, p oo, P coincides with (Tu) for D < 4. 
In the cut-ofF regularization the cut-off scale is used as 
the upper limit (instead of oo) in the second term of the 
right hand side in Eq. ([55]) . Since the last term has non- 
vanishing value, P does not coincide with (Tu). 



V. RADIUS AND MASS OF DENSE STAR 

To see the physical implication of the regularization de- 
pendence we evaluate the radius and the mass of dense 
star in both the dimensional and the cut-off regulariza- 
tions. For this purpose we confine ourselves to the case 
of a static and spherically symmetric star and do not 
care about the contribution of the strange quark and the 
neutrality of the star. 

Here we numerically evaluate the TOV equations in 
four dimensions using the results of the previous sections. 
Inside the stars the gravitational force should balance the 
pressure of the matter. For a static and spherically sym- 
metric star this condition is expressed by TOV equation, 



dPjr) 
dr 



dM{r) 
dr 



-G- 



p{r) + P{r) 



r{r - 2GM{r)} 
x{Af(r) +47rr^P(r)}, 

.2, 



= 47rr /9(r), 



(39) 
(40) 



where r is the radial distance from the center of a star. 
The pressure and the energy density are defined by 
the energy-momentum tensor P{r) — {Tu) and p{r) — 
— (T44). Because of the spherical symmetry the pressure, 
the energy density and the mass function M{r) are given 
as functions of r. 

We assume that the color superconducting phase can 
be realized inside dense stars and solve the Eqs. and 
(110)) with account of the EoS based on the extended NJL 
model [131 ■ We set the initial condition at r = as 
M(r = 0) = and p{r = 0) ^ and integrate the TOV 
equations. The energy density and the pressure decrease 
monotonically from the center to the surface. In our nu- 
merical analysis we define the surface of the quark star at 
the energy density, p ~ 0.0017GeV^, which corresponds 
to the ordinary energy density inside protons. 

In Fig. 10 we plot the mass of a star as a function of 
the central energy density. For lower pc the mass, M, 
increases with the central energy density. Increasing pc, 
we observe the maximal value for M. The mass decreases 
after passing through the maximal values. No stable ob- 
ject can be constructed in this parameter range with a 
negative slope, as is shown by dotted lines in Fig. 10. 
The maximal value for the mass in the dimensional reg- 
ularization is a few times as heavy as that in the cut-off 
regularization. 

In Fig. 11 the mass of the star is plotted as a function of 
the radius r. The dotted lines correspond to the unstable 
solution in Fig. 10. In the dimensional regularization we 
obtain a heavier and a larger quark star. It means that 



1.5 



1.0 



0.5 



1 ' 

dimensional 


1 1 1 1 1 M| 1 .1 


regularization 


- 




cut-off 


,,/,,, ,1 


regularization 





0.0001 



0.001 0.01 
Pc (GeV) 



FIG. 10: Mass of the star as a function of the energy density 
at the center of the star. Solid and dotted lines denote the 
stable and unstable states respectively. 
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FIG. 11: Relationship between mass and radius of the star. 
Solid and dotted lines denote the stable and unstable states 
respectively. 



the solution of TOV equations has a strong dependence 
on the regularization scheme. Since we do not include the 
contribution from the strange quark, our result obtained 
in the cut-off regularization shows a qualitatively similar 
behavior but does not exactly reproduce the result of the 
Ref. M. 



VI. CONCLUDING REMARKS 

We have investigated the phase structure of the chi- 
ral symmetry breaking and the diquark condensation 
(5J7c(3) symmetry spontaneous breaking) on the basis 
of the extended two-flavor NJL model by using the di- 
mensional and the cut-off regularizations. 

Evaluating the effective potential, we obtain the ex- 
pectation values of a and A. The values show a different 
behavior for the different regularizations. In the dimen- 
sional regularization we obtain a larger mass gap and no 
coexistence phase with non-vanishing (a) and (A) at the 



9 



quark massless limit, m ^ 0. The diquark condensation, 
(A) 7^ 0, is caused by the Cooper instabiUty of quark 
pairs near the Fermi surface. We loose such a contribu- 
tion in the cut-off regularization, as the chemical poten- 
tial approaches to the cut-off scale. The diquark con- 
densation is suppressed for a larger chemical potential, 
fi A. In the dimensional regularization |(A)| monoton- 
ically increases as the chemical potential increases. 

Next we studied the phase structure and the quark 
number susceptibility. TCP and CEP for the dimen- 
sional regularization are observed at higher temperature 
and lower chemical potential than those for the cut-off 
regularization. As is seen in Fig. 5, a complex struc- 
ture is realized around CEP for the coexistence phase. It 
is well-known that in the framework of the perturbative 
QCD the 2SC is predicted even at asymptotically high 
values of chemical potential /i. However, one of the de- 
ficiencies of the NJL model with cutoff regularization is 
the absence of the 2SC phase at sufhciently high values 
of /i. In contrast, it is shown that in the NJL model with 
dimensional regularization the 2SC phase is presented in 
the phase diagram at rather large values of /li, and this 
conclusion is in agreement with perturbative QCD anal- 
ysis. 

The NJL model used in the paper does not describe 
confinement. For the future study of the confinement- 
deconfinement transitions we are aware of the papers 
where confinement is simulated in NJL by introduction 
of the IR cutoff (to lift up the threshold of the decay into 
quarks) ,50] . There is also another way to address 
confinement [5l|, where they just compare the binding 
energies per a quark for the case of diquark and nucleon 
clustering. 

Confinement may have been relevant to our previous 
paper where we discussed the width of sigma in 

the leading order. In this paper we claimed that 
the problem of the unphysical width can be fixed at the 
next order, where the decay of sigma into 2tt takes 
place. So we hoped to fix the problem of the sigma width 
without even simulating confinement. 

We have also evaluated the energy-momentum tensor 
and the EoS. The dependence on the regularization ob- 
served in the high density region is not small. As is seen 



in Fig. 9, the rate of (Tu) and — (T44) is near unity in 
the dimensional regularization. It has been shown that 
thermodynamic pressure does not coincide with the space 
component of the stress tensor (Tu) in the cut-off regu- 
larization. 

Evaluating the TOV equations, we have obtained the 
relationship between the mass and the radius of dense 
stars in both regularizations. The existence of heavier 
and larger quark stars is expected in the dimensional 
regularization. The regularization dependence is mainly 
generated by the difference of the mass gap in Figs. 1 and 
2. The diquark condensation causes only a little contri- 
bution for the mass and radius of dense stars [28|. It 
should be noted that the size of the mass gap depends 
on the regularization parameter D. We obtain a smaller 
mass gap for larger dimensions [37| . However, a lower 
dimensional four-fermion interaction model is favored in 
another approach to QCD phenomena. If we adopt the 
ladder and the instantaneous exchange approximations 
and neglect momentum dependent parts of the fermion 
self-energy, Schwinger-Dyson equation coincides with the 
gap equation in the two dimensional NJL model at the 
leading order of expansion (53 |. 

The present work is mainly restricted to the analysis 
of the regularization dependence of observables in the 
extended NJL model. To apply our result to a more 
realistic case of a quark star we should include the con- 
tribution of the strange quarks and take into account the 
electric and the color neutrality. It would be also inter- 
esting to study the masses of tt and a mesons in the 2SC 
phase using the dimensional regularization and to com- 
pare the results with those in the 2SC phase obtained in 
the extended NJL model with cut-off regularization [5^ . 
Currently we are working on these issues. 
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